Analytical approximation schemes for solving exact renormalization 
group equations in the local potential approximation 



C. Bervillier, B. Boisseau, H. Giacomini 

Laboratoire de Mathematiques et Physique Theorique, 
UMR 6083 (CNRS), 
Federation Denis Poisson, 
Universite Francois Rabelais, 
Pare de Grandmont, 37200 Tours, France 



Abstract 

The relation between the Wilson-Polchinski and the Litim optimized ERGEs in the local potential approximation is 
studied with high accuracy using two different analytical approaches based on a field expansion: a recently proposed 
genuine analytical approximation scheme to two-point boundary value problems of ordinary differential equations, 
and a new one based on approximating the solution by generalized hypergeometric functions. A comparison with the 
numerical results obtained with the shooting method is made. A similar accuracy is reached in each case. Both two 
methods appear to be more efficient than the usual field expansions frequently used in the current studies of ERGEs 
(in particular for the Wilson-Polchinski case in the study of which they fail). 
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1. Introduction 

The non-decoupling of the relevant scales on a wide and continuous range of magnitudes in many areas 
of physics has led to the invention (discovery) of the renormalisation group (RG) pQ. Whereas they have 
been discovered in the framework of the perturbative (quantum held) theory, the RG techniques tackle a 
nonperturbative physical phenomenon [SJ. Nonperturbative approaches are difficult to implement and to 
control, and during a long time one has essentially carried on perturbative RG techniques (see, e.g., [5]). 
Nowadays, the huge growth of the computing capacity has greatly modified this behaviour pattern and, 
already since the beginning of the ninety's, one has considered [J] with a greater acuteness the exact RG 
equations (ERGEs) originally introduced by Wilson [5], Wegner and Houghton [6] in the seventy's and 
slightly reformulated by Polchinski [7] in the eighty's (for some reviews on the ERGEs see [5]). 

Initially, the ERGEs are integro-differential equations for the running action S [(f), t] [assuming that (j> (x) 
generically stands for some held with as many indices as necessary and t — — ln(A/Ao) the logarithm of 
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a running momentum scale A]. They have been extended to the running (average) effective action T [<p,t] 
[914] . Such general equations cannot be studied without the recourse to approximations or truncations. One 
of the most promising approximations is a systematic expansion in powers of the derivative of the field 
(derivative expansion) [10] which yields a set of coupled nonlinear partial differential equations the number 
of which grows quickly with the order of the expansion. In the simplest cases (e.g., for the scalar field), 
the determination of fixed points (and of their stability) amounts to study ordinary differential equations 
(ODEs) with a two-point boundary value problem that may be carried out numerically via a shooting (or a 
relaxation) method. 

A pure numerical study is in general not easy to implement and to control. For example, in the shooting 
method, the discovery of the right adjustment of the parameters at the boundaries requires a good knowledge 
a priori of their orders of magnitude (initial guesses). It is thus interesting to develop concurrently some 
substitute analytical methods. A popular substitute to the ODEs of the derivative expansion is provided by 
an additionnal expansion in powers of the field which yields a set of coupled algebraic equations which may 
be solved analytically, at least with the help of a symbolic computation software. Various field expansions 
have been implemented with more or less success |ll|12|13|14j . Unfortunately, the methods proposed up to 
now, if they are easy to implement, do not work in all cases and especially in the most famous and simplest 
case of the Wilson-Polchinski ERGE [5|7] (equation for the running action S [<f>,t] with a smooth cutoff). 

The object of this paper is to present two new substitute analytical methods for studying ODEs which, at 
least in the local potential approximation of the derivative expansion (LPA), works for the Wilson-Polchinski 
ERGE. One of the methods, recently proposed in [15j, is a genuine analytical approximation scheme to two- 
point boundary value problems of ODEs. The other method is new. It is based on approximations of the 
solution looked for by generalized hypergeometric functions. It has a certain similarity with another new 
and interesting method based on the representation of the solution by Pade approximants just proposed in 
[16j by P. Amore and F. M. Fernandez independantly from the present work. We illustrate the effectiveness 
of the two methods with the explicit consideration of two ERGEs in the local potential approximation: 
the Wilson-Polchinski equation and the Litim optimized RG equation [17] for the running effective action 
(named the Litim equation in the following). Following a conjecture first stated in |18|19j . the equivalence of 
these two equations (in the LPA) has been proven by Morris [5D] and recently been numerically illustrated 
[21j with an unprecedented accuracy for the scalar field in three dimensions (d — 3). This particular situation 
provides us with the opportunity of testing efficiently the various methods of study at hand. 

The following of the paper is divided in five sections. In section [21 we briefly present the direct numerical 
integration of the ODEs for the scalar model using the shooting method: determinations of the fixed point 
and the critical exponents for both the Wilson-Polchinski and Litim equations in the LPA (distinguishing 
between the even and odd symmetries) . A brief presentation of the currently used field expansion is given in 
section [3] In section [4] we analyse several aspects of the method of [15] applying it to the study of the two 
equations. We calculate this way the fixed point locations with high precision and compare the results with 
the estimates obtained in section [2] We show how the leading and the subleading critical exponents may 
be estimated using this recent method. In section [5] we present a new approximate analytical method for 
ODEs which is based on the definition of the generalized hypergeometric functions. We show that it is well 
adapted to treat the Wilson-Polchinski case whereas the Litim case is less easily treated. We relate these 
effects to the convergence properties of the series in powers of the field. Finally we summarize this work and 
conclude in section [6] 



2. Two-point boundary value problem in the LPA 

In this section we briefly present the two-point boundary value problem to be solved in the LPA of the 
ERGE. The Wilson-Polchinski equation is first chosen as a paradigm in section |2~T1 The principal numerical 
results obtained from the numerical integration of the ODE using the shooting method are given. In section 
(|2.2p . the Litim equation is also studied. 
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2.1. Wilson-Polchinski's flow equation for the scalar-field 



The original Wilson-Polchinski ERGE in the LPA expresses the evolution of the potential U (4>, t) as 
varying the logarithm of the momentum scale of reference t = — In (A/Ao) (with <p G M). In three dimensions, 
it reads: 

U = U" - (U'f - ~(f>U' + 3U , (1) 
in which U = dU (fa t) /dt, U' = dU (<p, t) /d<f>, U" = d 2 U (0, t) /d<j> 2 . 

2.1.1. Fixed point equation 

The fixed point equation corresponds to U = 0. It is a second order ODE for the function U (<f>): 

U" - (U'f - ^4>U' + 3U = , (2) 

the solution of which (denoted U* (cf>) below) depends on two integration constants which are fixed by two 
conditions. The first one comes from a property of symmetry assumed to b^3 U* (—</>) = U* (4>) which 
provides the following condition at the origin for U* ((f)): 

£T'(0) = 0. (3) 

The second condition is the requirement that the solution we are interested in must be non singular in the 
entire range <j> € [0, oof. Actually, the general solution of ([2} involves a moving singularity [22] of the form: 

losing = - In 1 O - cf>\ , (4) 

depending on the arbitrary constant <j> . Pushing <p Q to infinity allows to get a non-singular potential since, in 
addition to the two trivial fixed points U* = (Gaussian fixed point) and U* = — f + 2 ( n i§ n temperature 
fixed point), eq. admits a non-singular solution which, for <fi — > 00, has the form: 

„ ... (j? L e I8b 2 ^ 1 108 b 3 4/5 \ 

U asy W = ^r+b^ + V -- + t + OU 4/5 , (5) 

z zo o 625 4> 5 v ' 

in which b is the only remaining arbitrary integration constant. The non trivial (Wilson-Fisher [23j ) fixed 
point solution which we are interested in must interpolate between eqs. ([3]) and ([5]). Imposing these conditions 
fixes uniquely the value b* of b which corresponds to the fixed point solution we are looking for. 

We have determined b* by using the shooting method }24| : starting from a value <p a supposed to be large 
where the condition (0 is imposed (with a guess, or trying, value of b ~ 6*), we integrate the differential 
equation ([2]) toward the origin where the condition ([3]) is checked (shooting to the origin), we adjust the 
value of b to b* so as the latter condition is satisfied with a required accuracy. A study of the stability of 
the estimate of b* so obtained on varying the value 4> a provides some information on the accuracy of the 
calculation. 

Rather than ([5]), it is more usual to characterize the fixed point solution from its small field behaviour: 

UW _ » . »V + Ml+jjO^, _ * (1 + 3*H1 ± 24*) ^ e + Q w (6) 

and to provide the value of either of the two (related) quantities: 

k* = U* (0) , (7) 
r* = U*" (0) = -3fc* . (8) 

In the shooting-to-origin method, the determination of r* (or k*) is a byproduct of the adjustment of b*. 
The adjustment of b* may be bypassed by shooting from the origin toward (j> a , then r* is adjusted in such 
a way as to reach the largest possible value of <p a . In that case b* is a byproduct of the adjustment. 



The other possibility U* (—(f) = —U* (<fi) gives only singular solutions at finite <j>. 
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-0.228 598 202 437 022 -2.296 3 10 
-0.228 598 202 4370219 -2.3116 20 
-0.228 598 202 4370219 -2.316 2 40 

Table 1 

The fixed point parameter r* is already well determined for rather small values of cj> whereas b* [fixed point value of b in JSJ] 
still is not. 

Because the boundary condition at <p a is under control, the shooting-to-origin method provides a better 
determination of r* than the shooting-/rom-origin method. However, this latter method is more flexible and 
may easily yield a rough estimate on r* which can be used as a guess in a more demanding management of 
the method. Notice that, due to the increase of the number of adjustable parameters, this way of determining 
a guess is no longer possible in a study involving several coupled EDOs. Consequently, the development of 
other methods as, for example, those two presented below is useful to this purpose (see also [16]). 

Table Q] displays the determinations of r* and b* for three values of 4> a - O ne ma y observe that a high 
accuracy on r* is required to reach a yet small value of <p a whereas b* is only poorly determined. Obviously, 
considering higher values of <j) a and/or higher order terms in eq. ((5|) allows to better determine b* , one more 
term in ([5]) and <f) a = 1000 yields: 

b* = -2.318 29, (9) 
but the estimate of r* is not improved compared to the values given in table [1] (the machine-precision was 
already reached). We finally extract from table [1] our best estimate of r* (or k*) as obtained from the study 
of the fixed point equation @ alone: 

r* =-0.228 598 202 437 022 ± 10 -15 , (10) 
k* = 0.076 199 400 812 340 7± 10~ 16 . (11) 

Individually, these values do not define the potential function U* (<fi) the knowledge of which requires the 
numerical integration explicitly performed in the shooting method. 



2.1.2. Eigenvalue equation 

The critical exponents are obtained by linearizing the flow equation |lj near the fixed point solution 
U* (</>). If one inserts: 

U(cf>,t) = U* (0)+ee At .g(0) , 
into the flow equation and keeps the linear terms in e, one obtains the eigenvalue equation: 

g"-2g'U*'-^g' + (3-X)g = 0. (12) 

Again it is a second order ODE the solutions of which are characterized by two integration constants. 

Since U* {<j>) is an even function of </>, eq. (fl~2|) is invariant under a parity change. Then one of the integration 
constants is fixed by looking for either an even or an odd eigenfunction g (<f) which implies either g' (0) = 
(even) or g (0) =0 (odd). The second integration constant is fixed at will due to the arbitrariness of the 
normalisation of an eigenfunction. Thus, assuming either g (0) = 1 (even) or g' (0) = 1 (odd), the solutions of 
(|12p depend only on A and on the fixed point parameter k* . For example, these solutions have the following 
expansions about the origin <f> = : 



#even {4>) 



(A -3) 

—2- ' 



(A -2 - 12 fc*) 



0\ 



12 

U\ - ■ (2A-5-12fc*) ^ 3 . 5 . 

When the fixed point solution U* is known, the values of A [the only remaining unknown parameter in (Q3 
are determined by looking for the solutions which interpolate between either g' (0) = (even) or g (0) = 
(odd) and the regular solution of (|12p which, for <\> — > oo, is: 
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V 


ft* 




0.649 561773 88011 


-2.318145 


12 


0.649 561 773 880 80 


-2.318 257 


22 


0.649 561 773 880 65 


-2.318 280 


32 


0.649 561773 880 65 


-2.318 285 
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Table 2 

Values of the critical exponent v determined together with ft* (and thus r*) whereas cf> a is varied. Compared to table [T] a 
better determination of 6* is obtained [see the best value of 6* given by eq. J9}] ■ 



0.655 745 939193 3.180 006 512 059 5.912230612 8.796092825 11.798087 66 14.896053176 

Table 3 

Best estimates of the six first subcritical exponents for the Ising-likc scalar model (i.e. even case), all digits are significant. 



9aey{<t>) = S </)~ 



1 + (3 - A) 



12b* 

25c0 



36b* 2 (2 A -3) 2 (2A-1) 



625. 



125 < 



O 



(13) 



in which b* is given by The value of So is related to the choice of the normalisation of the eigenfunction 
at the origin, it is a byproduct of the adjustment in a shooting-/rom-origin procedure. 

In the even case, it is known that the first nontrivial positive eigenvalue Ai (there is also the trivial value 
Ao = d = 3), is related to the critical exponent v which characterizes the Ising-like critical scaling of the 
correlation length £. One has v — 1/Ai and the first negative eigenvalue, A2, is minus the Ising-like first 
correction-to-scaling exponent u>x (u>x — — A2) and so on. 

In the odd case, the two first (positive) eigenvalues are trivial in the LPA. One has: 



Ai = 



A 2 = 



d + 2-r] 
2 

d-2 + T] 



(14) 
(15) 



in which 77 is the critical exponent which governs the large distance behaviour of the correlation functions 
right at the critical point, it vanishes in the LPA. With the dimension d — 3 and the approximation 
(LPA) presently considered, (fl4|) and (TlB"]) reduce to Ai = 2.5 and A2 = 0.5. Consequently the first non- 
trivial eigenvalue is negative and defines the subcritical exponent #5 = fox — —A3 sometimes considered to 
characterize the deviation of the critical behaviour of fluids from the pure Ising-like critical behaviour. 

To determine the eigenvalues we use again the shooting-to-origin method with the two equations (J2J [T2|) . 
However, in addition to A, we leave also b* adjustable instead of fixing it to the value given in ©. 

In the even case, the values we obtain for v and b* are shown in table [5] for four values of <j> a . Comparing 
with the values displayed in table [T] one observes a better convergence of b* to the best value (O whereas r* 
remains unchanged compared to (JTUJ) . As for the best estimate of v, it is: 



z, bcst = 0.649 561 773 880 ± 10" 



(16) 



that is to say: 

Aibcst = 1.539 499 459 808 ± 10~ 12 . (17) 

We have proceeded similarly to determine the Ising-like subcritical exponent values displayed in table [31 
In the odd case, we obtain: 

u)x = 1.886 703 838 091 ± 10~ 12 . (18) 

Table [4] displays the values of the other subcritical exponents of the same family as u> but with a lower 
accuracy. Of course, the values presently obtained are in agreement with the previous estimates |25I21] . 
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U)2 u>4 

4.524 390 734 7.337650 643 10.283 900 73 

Table 4 

Best estimates of the odd-case subcritical exponents other than ui\ for the scalar model. 

2.2. Litim's flow equation for the scalar field 

Following a conjecture first stated in [18ll9j . the equivalence in the LPA between the Wilson-Polchinski 
flow ((1]) and the Litim optimized ERGE [17] for the running effective action T [tp, t] has been proven by 
Morris [20]. The Litim flow equation for the potential V (ip,t) reads in three dimensions (compared to [20] 
an unimportant shift V — ► V — 1/3 is performed): 

v = 1 -TTv^-^ v ' + 3V - (19) 

It is related to (JTJ) via the following Legendre transformation: 

1 



^4> 2 -u(<t>,t) 

<p = <t>-U'{ct>,t) 



(20) 



The general solution of the fixed point equation (V = 0) involves the following moving "singularity" (V" 
is singular) at the arbitrary point <p : 

%i„ g (^ = -i + T 4=|^ -¥'| 3/2 ■ (21) 

2.2.1. Fixed point solution 

The numerical study of the fixed point solution of (fl"9)) follows the lines described in the preceding sec- 
tions. This may be done independently, but due to ([20]) . one may already deduce from the previous study 
the expected results. Similarly to ([5]), the asymptotic behaviour of the non trivial fixed point potential is 
characterized by the integration constant 6l in the following expression [deduced from IT9|) ]: 

= b L ^ - \ + 1 4 - 1 + Q (y- 12 ) . (22) 
3 150 ol</5 4 6300 

It is easy to show from and ([2U)l that the value b* L we are looking for is related to b* as follows: 

b* --1(1 " 
then, from the previous result © we get: 

&2 - 0.001000 25. (23) 

Similarly for the potential parameters 

k* L = V* (0) , 
r* L = V*" (0) , 

which correspond to b* L , they are related to the Wilson-Polchinski counterparts k* and r* as follows: 

k* L =k*, (24) 

r * L =r^- (25) 

This latter relation, using ([T0|) . gives: 

r£ ~ -0.186 064 249 470 314 ± 10~ 15 . (26) 
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As precedingly, those values do not provide the potential function V* (if) the knowledge of which requires 
an explicit numerical integration. 



2.2.2. Eigenvalue equation 

A linearization of the flow equation (fT9|) near the fixed point solution V* (<p): 

V{y,t)=V* (<p)+ee xt h(<p) , 



provides the Litim eigenvalue equation: 



(3 - A) h — + 

2 (l + V* 



= 0. 



Taking into account 
form: 



one can show that (|2~T|) admits a regular solution which, for if 

1 1 



h asy (if) = W (3 ~ A) < 1 - (A — 3) (2 A — 5) 



■o(<p- 



(27) 
oo, has the 

(28) 



2250 bfip w 47250 bf ip u 

in which b* L is given by (|23[) . In the following we may set Si = 1 since the normalisation of the eigcnfunction 
may be chosen at will. 

As precedingly, we must distinguish between the odd and even eigenfunction h (if). The shooting method 
gives the same values as in the Wilson- Polchinski case (see |18|26|27|21j ) and we do not present them again. 



3. Expansion in powers of the field 



In advanced studies of the derivative expansion 28J or other efficient approximations of the ERGE [29] and 
in the consideration of complex systems via the ERGEs |30j . a supplementary truncation in powers of the 
field is currently used (see also [8]). With a scalar field, this expansion transforms the partial differential flow 
equations into ODEs whereas the fixed point or eigenvalue ODEs are transformed into algebraic equations. 
Provided auxiliary conditions are chosen, the latter equations are easy to solve analytically using a symbolic 
computation software. Actually the auxiliary conditions currently chosen are extremely simple: they consist 
in setting equal to zero the highest terms of the expansion so as to get a balanced system of equations. 

A first kind of expansion, about the zero field -referred to as the expansion I in the following, has been 
proposed by Margaritis et al [TT] and applied to the LPA of Wegner-Houghton's ERGE [5] (the hard cutoff 
version of the Wilson-Polchinski equation). A second kind of expansion, relative to the (running) minimum 
of the potential (expansion II), has been proposed by Tetradis and Wetterich [12j and more particularly 
presented by Alford [13] using it, again, with the sharp cutoff version of the ERGE. 

It is known that, for the Wegner-Houghton equation in the LPA, expansion I does not converge due to 
the presence of singularities in the complex plane of the expansion variable [31] . Expansions I and II have 
been more concretely studied and compared to each other by Aoki et al in [14] who also propose a variant 
to II (expansion III) by letting the expansion point adjustable. They showed, again on the LPA of the 
Wegner-Houghton equation, that expansion II is much more efficient than expansion I although it finally 
does not converge and expansion III is the most efficient one. Expansions II and III work well also on the 
ERGE expressed on the running effective action (effective average action, see the review by Berges et al 
in [5]). The convergence of those expansions have also been studied in [5S] according to the regularisation 
scheme chosen and in particular for the Litim equation (|19p . In this latter study it is concluded that both 
expansions I and II seem to converge although II converges faster than I. 

A striking fact emerges from those studies, the Wilson-Polchinski equation in the LPA, the simplest 
equation, is never studied using the field expansion method. The reason is simple: none of the expansions 
currently used works in that case. 

Actually the strategy of these methods, which consists in arbitrarily setting equal to zero one coefficient 
for the expansion I and two for the expansions II and III, is probably too simple. With regards to this kind 
of auxiliary conditions, the failure observed with the Wilson-Polchinski equation is not surprising and, most 
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certainly, there should be many other circumstances where such simple auxiliary conditions would not solve 
correctly the derivative expansion of an ERGE. 

In the following sections we examine two alternative methods with more sophisticated auxiliary conditions. 
We show that they yield the correct solution for the Wilson-Polchinski and its Legendre transformed (Litim) 
equations. Both methods are associated to expansion I (about the zero- field). The first one has recently been 
proposed in [15j as a method to treat the two point boundary value problem of ODEs. It relies upon an 
efficient account for the large field behaviour of the solution looked for. An attempt of accounting for this kind 
of behaviour within the field expansion had already been done by Tetradis and Wetterich via their eq. (7.11) 
of |12j . In the present work, a much more sophisticated procedure is used. It relies upon the construction of an 
added auxiliary differential equation (ADE). We refer to it in the following as the ADE method. The second 
method is new. It relics upon the approximation of the solution looked for by a generalized hypergeometric 
function. We refer to it in the following as the hypergeometric function approximation (HFA) method. 

4. Auxiliary differential equation method 

Let us first illustrate the auxiliary differential equation (ADE) method on the search for the non trivial 
fixed point in the LPA for both the Wilson-Polchinski equation ([2]) and the Litim optimized equation (|T9|) . 
Since there are two boundaries (the origin and the "point at" infinity), we distinguish between two strategies. 

- An expansion about the origin in the equations (small field expansion) and the account for the leading 
high field behaviour of the regular solution which we are looking for. This determines the value of r* or 

- A change of variable (f> — ► l/(f> or ip — ► 1/ip which reverses the problem: an expansion about infinity (new 
origin) in the equations (high field expansion) and the account for the leading small field behaviour of the 
regular solution which we are looking for. This determines the value of b* or b* L . 

4.1. Wilson-Polchinski's fixed point 

4.1.1. Small field expansion and leading high field behaviour 

For practical and custom reason^], instead of ((2]) we consider the equation satisfied by the function w (x) 
related to the derivative of the potential U' (<fi) as follows: 

U'(cj>) =<M0 2 ) . (29) 
so that, with x = 4> 2 , the fixed point equation ||3J| reads: 

4 x w" - 2 w 2 - 4 x w w' + (6 - x) w' + 2 w = , (30) 

in which a prime indicates a derivative with respect to x. 

This second order ODE has a singular point at the origin and, by analyticity requirement, the solution 
we are looking for depends on a single unknown integration-constant (noted r below). 

Let us first introduce the expansion I of Margaritis et al |llj . The function w (x) is expanded up to order 
M in powers of x: 

M 

wm {x) = r + y"] a n x n , (31) 
n=l 

and inserted into the fixed point equation (|30p . 

Requiring that (|30p be satisfied order by order in powers of x provides an unbalanced system of M 
algebraic equations with M + 1 unknown quantities {r, ai,- ■ ■ ,om} [eq. ([30|) is then satisfied up to order 
M — 1 in powers of x] . With a view to balancing the system, om = is simply set equal to zero and if the 
solution involves a stable value as M grows, then it constitutes the estimate at order M of the fixed 

2 The change x = (j> 2 is useful in practice to avoid some degeneracies observed in 1 151 when forming the auxiliary differential 
equation. Taking the derivative / = U' is only a question of habit. 
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point location corresponding to expansion I. As already mentioned, in the case of the Wilson-Polchinski 
equation (|30|) under study, the method fails: all the values obtained for r* RI are positive whatever the value 
of M whereas the correct value should be negative as shown in section [2. 1.1 1 

In the ADE method, the condition a a/ = is not imposed. The previous algebraic system is first solved 
in terms of the unknown parameter r so as to get the generic solution of (|30p at order M in powers of x: 

M 

w M (r;x)=r + ^2 a n (r) x n . (32) 

n=l 

In order to get a definite value for r, instead of arbitrarily imposing om (r) = 0, an auxiliary condition 
is formed which explicitly accounts for the behaviour at large 4> given by ([5]) . With w (x) , this behaviour 
corresponds to: 

Wasyix) = 1, (33) 
x— »oo 

w' asy (x) = 0. (34) 

The auxiliary condition is obtained via the introduction of an auxiliary differential equation: 

- Consider a first order differential equation for w (x) constructed as a polynomial of degree s (eventually 
incomplete) in powers of the pair (w, w'): 

Gi +G 2 w + G 3 w' + G 4 w 2 + G 5 w w' + G 6 w' 2 +••• + £„ w s - q w' q = , (35) 

in which, when the degree s of the polynomial is saturated then q = s and the number n of coefficients 
Gi is equal to (s + 1) (s + 2) /2, conversely when it is not then < q < s and n = s(s + l)/2 + q + 1. 

- The constant coefficients Gi are then determined as functions of r by imposing that the solution wm (r; x) 
of (|30|) previously determined for arbitrary r at order M in powers of x be also solution of (|35|) (at the 
same order M). Due to an arbitrary normalisation which allows to fix, for example G\ = 1, a simple 
counting shows that the identification implies M = n — 1. The resulting set {Gi(r);i — 2,...,n} is 
formed of rational functions of the unknown parameter r. Hence, a new differential equation for w (x) is 
obtained: 

1 + G 2 (r) w + G 3 (r) w' + G 4 (r) w 2 + G 5 (r) w w' + G 6 (r) w' 2 + ■ ■ ■ + G n (r) w g - q w' q = 0, (36) 

which is satisfied by construction at order M in powers of x by (I32p which is already solution at the same 
order of (|3"0j) . 



- The last step is then to impose that the new equation (|36p be also satisfied when x — > oo. Taking into 
account (f33l l34f it comes the final auxiliary condition: 

1 + G 2 (r) + G 4 (r) + ■ ■ ■ + G s(s+1)/2+1 (r) = . (37) 

Solving this auxiliary condition for r amounts to determining the roots of a polynomial in r. As the order 
M grows some root values appear to be stable. Those stable values are candidates for the fixed point solutions 
we are looking for. In a way similar to |16j . the obtention of the auxiliary condition may be obtained without 
determinating explicitly the coefficient functions Gi (r). For this, it is sufficient to consider the matrix T 
of the homogeneous system of linear equations for all the G^s formed with eq |35|) to which is added its 
expression when x — * oo. When the function w (x) is replaced by the expansion (|32|) at the required order 
the matrix T depends only on the coefficients a n (r) of the Taylor expansion (132|) and the auxiliary condition 
then finally reduces to: 

detT = (38) 

Before going further, it is worthwhile indicating that a variant of the method which consists in remplacing 
w' by xw' in the auxiliary differential equation (|35[) has appeared more efficient [e.g., see figure [5]). 

Figure ([T]) shows the distribution of all the real roots Tm of (|3"T|) for the variant as the order M varies up 
to 28. The three expected fixed points encountered in section (|2 . 1 . 1[) are clearly evidenced by a threefold 
accumulation about the respective values 1 (HT), (Gaussian) and r* (Wilson-Fisher). Although a huge 
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Fig. 1. Distribution of the real roots tm (open circles) of (1376 as function of the order M of the Taylor series about the origin 
d 3 1 1 ) [with the ADE pair (w,xw')]. A threefold accumulation occurs about the expected fixed points: trivial high temperature 
(r* = 1) and Gaussian (r* = 0) fixed points and about the non-trivial Wilson-Fisher fixed point (r* ~ —0.2286) [LPA, d=3]. 



accumulation of roots around the right value occurs, the approach to r*, which we are interested in, may be 
followed step by step as the order M grows. 



4.1.1.1. Selection of the root To select the right value of the root corresponding to the nontrivial Wilson- 
Fisher fixed point, the following procedure has been applied. We know that the root of interest is negative 
and real, then we select the first negative real root that appears at the smallest possible order. At the next 
order we choose the real root the closest to the previous choice and so on. We obtain this way with M = 28 
the following excellent estimate: 

r* = -0.228 598 202 437 02 , (39) 

which coincides, up to the 14" 1 digit, with the estimate l|10p obtained by the shooting method. Figure [2] 
shows the accuracy obtained on r* by selecting the roots this way as M varies. 
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Fig. 2. Approximate number of accurate digits = — log |1 — t'm I obtained on the selected roots rjwr as a function of M 
and for two ADE pairs; the original (w,w'): crosses, and the variant (w, xw'): open circles [r* is given in eq. 110H . The highest 
values of M in each case is limited by time computing. A better efficiency is obtained with the variant. 

4.1.2. Subleading high field behaviour 

Equations (|33l [34)) used in the preceding calculations express exclusively the limit of the solution w* (x) 
when x — » oo, and we get the unique condition (|37p to estimate r*. In fact there are higher correction terms 
to (|331 134[) which vanish as x — > oo [the first of which correspond to those written in ([5])] . Such subleading 
contributions may as well be imposed in (|36p . In so doing, we require the auxiliary differential equation to 
be satisfied not only at infinity but also in approaching this point. Consequently we obtain several auxiliary 
conditions similar to (|37p . each of them corresponding to the cancellation of the coefficient of a given power 
of x. We have used them to determine r* again (the asymptotic constant b factorizes in the first subleading 
conditions so obtained). The results are similar to those obtained precedently with the leading conditions 
(f33"l [34]) alone. We have observed only a slight decrease in the accuracy: the higher the subleading term 
considered the weaker the convergence to r* . This shows the coherence of the ADE method: the auxiliary 
condition is not an isolated point condition, it emanates from a differential equation constructed to be 
satisfied by the function looked for. 

When the order of the subleading contribution is high enough, the constant b no longer factorizes and the 
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subleading auxiliary condition depends non trivialy on the (non-independent) integration constants (r and 
b) characterizing the fixed point solution. We have tried to determine the value b* by imposing the individual 
vanishing of such contributions for r — r* . Unfortunately, at the orders considered, the only knowledge of r* 
suffices to satisfy the condition (whatever the value of b) . It is possible that considering much higher orders 
would allow us to get an estimate of b* this way. 

4.1.3. High field expansion and leading small field behaviour 

With the determination of b* by the ADE method in view, let us perform the change of variable x — > y~ 5 
and the following change of function: 

u(y)=y- 2 [w(y- 5 )-l] , (40) 
so that, from (0 and ([29]) . u* (y) has the following form for small y : 

u* (y) = A* + ±A* 2 y 2 - 1a*V + O (y 5 ) , (41) 

with 

A* = 66*/5 . (42) 
The fixed point differential equation ^ is then transformed into: 

- 10 yu 2 + 5 (5 + 2y 5 ) u' - 4 (y 4 -5y 2 u') u + Ay e u" = 0, (43) 
the solution of which must satisfy the following condition, see (|4ip : 

u* (0) = A* , 
u*'(Q) = 0, 

with A* to be determined so as, using ([3T1 [40]) . to get at infinity: 

< sy (y) = o . 

The ADE method described in the preceding sections is used to determine the value of A*. Since there 
are some holes in the first terms of the series (|41[) . the first significant estimates are obtained for values of 
M higher than in section 14.1.11 Figure shows that the selected sequence of roots corresponding to A* 
converges to —2.73532 whereas, according to (f9l [42]) , the right value expected from the shooting method is 
-2.78195. 

This failure of the ADE method in determining correctly A* is presumably due to the zero radius of 
convergence of the Taylor series of u* (y) about y = 0. Actually, we have estimated this radius as the limit of 
the ratio of two consecutive terms and observed that it goes slowly but continuously to zero as the order M 
increases. This contrasts with the case of w (x) for which the same procedure quickly tends to the following 
finite limit for the fixed point solution corresponding to (fTOl [39]) : 

Rwp = 5.72167. (44) 

Notice that, although the ADE method does not provide the right estimate of A* (or b*), it gives a value 
close enough to it to be used as a guess in the shooting method. 

4.2. Litim's fixed point 

4.2.1. Small field expansion and leading high field behaviour 

For convenience we perform the following change, compared to section [2. 21 

V(<p)=w(<p 2 )-±, (45) 

so that the fixed point equation corresponding to (|19p reads (with x — ip 2 ): 
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Fig. 3. In the Wilson-Polchinski case, A* does not converge (dots) to the right value —2.78195 (horizontal line) but to —2.73532. 



3w-xw' = 0. (46) 

1 + 2 w> + 4 x w" 

The singularity at x = of this second order ODE allows us to look for an analytic solution which satisfies, 
in terms of a single unknown parameter fc, the following conditions at the origin: 

w(Q) = k, (47) 

•'M-a-s- (48) 

with fc adjusted to fc* so as to reach at infinity [from (f22l) ]: 

^ - b ^ + 150^ " + ° ( ^ 6) • (49) 

The expected value of fc* is related to r* L given in (|26|) as: 

^ fc* = 1 



3(l+rJ) 
13 



It is also related to k* given in (TTT]) via ([Ml |4"5")) as fc* = k* + 1/3. Consequently the estimation by the 
shooting method is: 

k* = 0.409 532 734 145 674 ± 1CT 15 . (50) 

The object of this section is thus to test whether the ADE method yields that value of k* [and also that 
of b* L given in [g5J) ]. 

Contrary to the Wilson-Polchinski case, the asymptotic behaviour (|4"5|) does not reach a finite value when 
x — > oo. But the third derivative of w* does. Hence, since 6^ I s still supposed unknown, the auxiliary first 
order differential equation (|35[) may be used with w and w' replaced respectively by w)' 4 ' and (where 
w^ n > stands for d n w/dx n ). Actually, both of these two derivatives go to zero as x — > oo so that finally the 
auxiliary condition similar to (|37|) . but with another normalisation of the Gi's (e.g. Gi = 1), reduces to: 

G 1 (k)=0, (51) 

whereas the function w (x) is expanded up to order M in powers of x and inserted into (|46[) to get the 
solution at this order as function of k: 

M 

w M (k;x) =k + y~l a n (Jo) x n . (52) 
n=l 

Similarly to the Wilson-Polchinski case, the complete set of real roots of (|5ip shows accumulations about 
the expected fixed point values. However the selection process described previously fails in picking the right 
value k* (of the nontrivial fixed point) although it is present among the roots. Actually, for M = 14 the 
selection gives 0.409 627 819 729 71 whereas a better value (0.409 532 733 212 35) exists at the same order 
[compare with (|50p ]. The variant utilised in the preceding case which consists in replacing by iw' 5 ' 
does not circumvents this difficulty. 

If instead of (w^ 4 \ w^) as ADE pair, we consider the combination h — 3w — xw' and its derivative hi 
with respect to x, (or the variant xhl to save some time computing), then the new pair, according to (|49p . 
vanishes also as x — > oo , and we observe, this time, that the selection process works again. This way, at 
order M = 19 the selection gives: 

k* = 0.409 532 73416, 

a value which coincides with (|50|) up to the 10 th digit. No doubt that considering higher values of M would 
have improved the accuracy. We note that, as with Wilson-Polchinski's function, the radius of convergence 
of the Taylor series of w (x) about the origin is finite, and is about: 

R L ~11.5. (53) 

Let us specify however that, contrary to the Wilson-Polchinski case, the test of the ratio aj/aj+i of two 
consecutive terms of the Taylor series about the origin does not converge. We have obtained (|5"3"]) by explicitly 
performing a partial summation of the series and studying it as a function of x. Nevertheless, we have also 
observed that the ratio |ai/at+3| raised to the power 1/3, roughly converges to (l53|) . This remark will have 
some importance in section 15.41 

Since expansions I and II work in the Litim case (see [26]), we can compare the ADE method with those 
two methods. Figure ^ shows the respective accuracies obtained on k* with the three methods as functions 
of the order M of the field expansion. One sees that expansion II and the ADE method provide better 
results than expansion I (which likely does not converge) and that the ADE method is most efficient than 
expansion II (we have not studied expansion III). 

4.2.2. Subleading high field behaviour 

As in the case of Wilson-Polchinski's equation, the subleading terms in (|49ll may be used to impose the 
auxiliary condition not only at infinity but also in approaching this point whatever the value of x. We 
observe the same phenomenon as in section f4. 1.21 the higher the subleading term considered the weaker the 
convergence to k* whereas b* L cannot be determined by imposing the individual vanishing of the subleading 
contributions for k = k* . 
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Fig. 4. Approximate number of accurate digits Nd = — log 1 1 — /cm /k* | [with k* given by (150H as functions of M for the 
estimations of the Litim fixed point value of k using three methods: expansion I (black dots), expansion II (crosses), and ADE 
(open circles). A better efficiency is obtained with the ADE method. 

However, the fact that the asymptotic behaviour (|49[) is an integer power of x provides us with the 
oportunity of determining b* L from the knowledge of k* as a boundary limit (a point condition) . Actually, 
since — > &b* L when x —> oo, we may choose (w^ 3 \ as ADE pair [or the variant (w)^ 3 \ xw^)}, and 

for k fixed to k* solve for b* L the resulting auxiliary condition at infinity. The accuracy on b* L obtained this 
way is not as large as in the case of k*, nevertheless, for M = 31 we obtain the following estimation: 

6^- 0.001007, (54) 

which is rather close to the shooting value (f2"3"|) . We indicate also that rough estimates of b* L already suf- 
ficiently accurate to be used as guesses in the shooting method are obtained for small values of M, e.g.: 
0.000 989 for M = 11 or even 0.0012 for M = 5. 

4.2.3. High field expansion and leading small field behaviour 

With a view to determining b* L directly by the ADE method, we invert the boundaries by changing the 
variable x — > y _1 and by performing the following change of function: 
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u(y) =fw ( - ) , (55) 



so that, from (j49|) . we deduce that u* (y) has the following form for small y : 

r ®=^i4^ +0feT (56) 
The differential equation for u (y) is: 

y 4 + l8y {u'f-u' [30u + y 2 (l + 4u")] =0. 

The solution must satisfy the following condition at the origin y — [see (|56[) ] : 

u* (Q) = b* L , 
u" (0) = , 

with b* L to be determined so as, using |47l l48l [55|) . to get at infinity: 

ul sy (y) = k*f+(±;-^j f + 0(y). 

As previously, we use the ADE method with a view to determining the value of b* L . For this we consider, 
the pair (u^ 4 \ u^) which vanishes at infinity (y — ► oo). Since there are some holes in the first terms of the 
series about the origin, see (j56|) . the first significant estimates are obtained for values of M higher than with 
the original function w (x). Although the positive roots obtained for 6l (we know that b* L is positive) have 
the right order of magnitude compared to |23|) the apparent convergent sequences do not provide the right 
value. Again, as in the Wilson-Polchinski case, we think that the failure of the ADE method is due to the 
(observed) zero radius of convergence of the Taylor series for u (y) about the origin. 

4.3. Eigenvalue estimates 

Let us consider the eigenvalue problem with the ADE method. This time two coupled nonlinear ODEs have 
to be solved together (the fixed point equation and the linearisation of the flow in the vicinity of the fixed 
point). We can solve these two equations together as the order of the field expansion M grows or consider 
separately the eigenvalue equation after having solved the fixed point equation with some accuracy. With 
the aim to be short, we present only the latter possibility which illustrates well the property of convergence 
of the method. 



4.3.1. Wilson-Polchinski's eigenvalues 

4.3.1.1. Small field expansion and leading high field behaviour Using a change of eigenfunction, g — > v, 
similar to (|29[) for the fixed point function, it comes: 
in the even case: 

g'{cp) = ^{4> 2 ) , 

and eq. ((T3|) yields the following behaviour at large x — <j) 2 : 



in the odd case: 
and eq. (fT3]) gives: 



v asy (*) = 

g> (0)= v(0 2 ) 
Vasy (x) = ^^So* (1 - 2A)/1 ° 



1 + 



1 + O (^ 2/5 ) 
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Fig. 5. Accumulations of real roots (open circles) of the auxiliary condition H37H about eigenvalues as the order M of the series 
varies in the Wilson- Polchinski even case. From top to bottom: X± = 1/u (second horizontal line), A2 = —ui\ (third h. line) and 
A3 = —u>2 (fourth h. line). A simple criterion of choice allows to determine their estimates at M = 20, see the values in eqs. 
J57H590 . An accumulation also occurs about the spurious value 5.8 (first h. line). 

The arbitrariness of the global normalisation of the eigenfunctions allows to choose v (0) = 1 (even) and 
v' (0) = 1 (odd) corresponding respectively to some definite values of So. So defined, the functions v (x) and 
v' (x) vanish at infinity provided that A > —2 in the even case and A > 1/2 in the odd case. Hence one 
could expect that, with the simple condition at infinity: v = v' = imposed in the auxiliary differential 
equation, the ADE procedure will, at best, allow the determination of exclusively the leading (Ai = 1/v) 
and first subleading (A2 = — u>x) eigenvalues in the even case and of only the trivial eigenvalue Ai = —oj\ in 
the odd case [see the values of these quantities in eqs. ((T7l [T8t and tables (014])]. Actually it is better than 
that since, as M grows, we observe among the real roots of the auxiliary condition for A that a hierarchy of 
successive accumulations takes place about the right values of the leading and subsequent eigenvalues [see 
figure [5] . 

Within each of these accumulations of real roots, we have been able to follow without ambiguity a con- 
vergent sequence to the right estimate. At order M = 20 with the ADE pair (v, v') supposed to vanish at 
infinity, and r* fixed to the value given in we have obtained the following estimates in the even case 
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^ = 0.649 561773 86, (57) 
w x = 0.655 745 92, (58) 
uj 2 = 3.178, (59) 

where the number of digits has been truncated with regard to the accuracy of the estimates obtained [by 
comparison with (fT6|) and table [3] . We see that the accuracy decreases as the order of the eigenvalue grows 
but also that we obtain an estimate of 102 whereas for that value v does not vanish at infinity. 

The same kind of observations stands in the odd case. We take the opportunity to indicate that choosing 
the ADE pair (/, /') with / = ^~^ v-iv' makes / vanish for A > —3/2 and the procedure gives a better 
accuracy on #1 than with the pair (v,v'). This way we obtain the following estimate [at order M = 20, 
compare with ([18])] 

#1 = 1.886 718. 

We have also noted the presence of accumulations of real roots about spurious positive values of order 5.8 
in the even case and 3.77 in the odd case. 



4.3.2. Litim's eigenvalues 

The determination using the ADE method of the eigenvalues from the Litim flow equation follows the 
same lines as previously for the Wilson-Polchinski flow equation. We limit ourselves in this section to a brief 
presentation of the main differences encountered. 

4.3.2.1. Small field expansion and leading high field behaviour Compared to (|27|) . we perform a change of 
eigenfunction, h — > vx, according to the symmetry considered: 
- in the even case: 

h(<p) = v L ( v 2 ) , 
then eq. (|28|) yields the following behaviour at large x = ip 2 : 

VLasy(i) = S X i ^ [l+0(x~ 5 )] . 

in the odd case: 

h(tp) = <pv L (ip 2 ) , 

and eq. (f2"8"| gives: 

v Lasy (x) = S 1 xW 2 -V[l + 0(x- 5 )} . 

So defined, the two functions v^ (x) and v' L (x) vanish at infinity provided that A > 3 in the even case and 
A > 5/2 in the odd case (whereas the arbitrary global normalisation of the eigenfunctions allows to choose 
V£ (0) = 1 (even) and v' L (0) = 1 (odd) corresponding respectively to specific values of Si). 

Although it works, the original ADE pair (vl,v' l ) is not the most efficient choice to obtain estimates 
of the first nontrivial eigenvalues. A better choice appears to be the pairs (/ (x) , /' (x)) with / (x) = 
(3 — A) vl (x) — xv' L (x) in the even case and / (x) = (5/2 — A) v^ (x) — xv' L (x) in the odd case (they 
correspond to eigenfunctions which vanish as x — > 00 for more negative values of A). With these choices 
and k* = 0.409 532 734 145 7 we identify immediately the trivial eigenvalues Ao = 3 in the even case and 
Ai = 2.5, A2 = 0.5 in the odd case but also, for M = 20, we obtain good estimates of the nontrivial leading 
and first sublcading eigenvalues: 

v = 0.649 561 774, u)i = 0.655 745 5, w 2 = 3.180 008, uj 3 = 5.896, 

#1 = 1.886 703 7, #2 = 4.5241, 

where the numbers of digits have been limited with respect to the estimated accuracy [compare with (|16|). 
table [3] (even) and (fT8|) . table [4] (odd)]. For each eigenvalue, the successive estimates may be followed 
unambiguously step by step when M grows so that the right values may be easily selected following the 
rules defined preccdcntly. 

We notice also the presence of spurious convergences and especially in the even case to the value about 
5.8 already encountered with the Wilson-Polchinski case. 
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5. Approximating by hypergeometric functions (HFA) 



The ADE method is most certainly efficient in many cases but it is relatively heavy regarding the com- 
puting time whereas the current methods, when they work, are lighter. In addition, none of these methods 
provides a global solution to the ODE studied: they yield an approximate value of the integration constant 
but not a function as global approximation of the solution looked for. 

We propose in this section an alternative method which is lighter than the ADE method and which 
provides a global approximation of the solution of interest. This new method is based on the definition 
property of the generalized hypergeometric functions. Let us first review the definition and main properties 
of these functions. 



5.1. Generalized hypergeometric functions 



For x £ C, a series 5* = X)^°=o a nX n is hypergeometric (see for example [32]) if the ratio a n +i/a„ is a 
rational function of n, i.e. 

Qn+i _ P(n) 
a n Q (n) ' 

for some polynomials P (n) and Q (n). 

If we factorize the polynomials, we can write: 

ftn+i _ (n + Qi) (n + a 2 ) ■ • • (n + a p ) . . 

a n ^ a0 (n + (3 1 )(n + f3 2 )---(n + f3 q )(n + l)- 

The factor (n + 1) in the denominator may or may not result from the factorization. If not, we add it along 
with the compensating factor in the numerator. Usually, the global factor ao is set equal to 1. 

If the set {cti\ includes negative integers, then S degenerates into a polynomial in x. 

When it is not a polynomial, the series S converges absolutely for all x if p < q and for |x| < 1/ |oo| if 
p = q + 1. It diverges for all x ^ if p > q + 1 . 

The analytic continuation of the hypergeometric series S with a non-zero radius of convergence is called 
a generalized hypergeometric function and is noted: 

pF q (ax, - • • , a p ;/3 1( • • • , (3 ■ a x) = — S . 

p F q (x) is a solution of the following differential equation (for a a = 1): 

[6(6 + P 1 -l)'-'(B + q -l)-x(6 + ai)'--(e + ap)] P F q (x) = , (61) 

where 

9 = x± 
ax 

When p > 2 or q > 1, the differential equation (|6Tj) is of order max (p,q + 1) > 2. It is of second order 
when q = 1 and p — 0, 1 or 2. It is of first order when q = and p = 1 

2-Fi is currently named the hypergeometric function. A number of generalized hypergeometric functions 
have also special names: qFi is called confluent hypergeometric limit function and iF\ confluent hypergeo- 
metric function. 

In the cases p < q for fixed and {/3j}, p F q (x) is an entire function of x and has only one (essential) 
singular point at x — oo. 

For p — q + 1 and fixed {cti} and in non-polynomial cases, p F q (x) does not have pole nor essential 
singularity. It is a single- valued function on the x-plane cut along the interval [1, oo], i.e. it has two branch 
points at x = 1 and at x = oo. 

Considered as a function of {Pf, i = 1, • • • , g}, p F q (x) has an infinite set of singular points: 

(i) (3 i — —m, m € N which are simple poles 

(ii) (3 i — oo which is an essential singular point (the point of accumulation of the poles). 
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As a function of {a^; i = 1, • • • pi 7 ^ (a;) has one essential singularity at each cxi — oo. 

The elementary functions and several other important functions in mathematics and physics are expressible 
in terms of hypergeometric functions (for more detail see [32 ) . 

The wide spread of this family of functions suggests trying to represent the solution of the ODEs presently 
of interest in this article, under the form of a generalized hypergeometric function. 



5.2. The HFA method 



For the sake of the introduction of the new method, let us first consider the Wilson-Polchinski fixed point 
equation (|30|) and the truncated expansion (|32f in which the coefficients a n (r) (n = 1, • • • , M) are already 
determined as function of r via a generic solution of (|3H)) truncated at order M (in powers of x). The question 
is again to construct an auxiliary condition to be imposed with a view to determining the fixed point value 
r* . To this end, by analogy with the generalized hypergeometric property definition recalled in section 15. li 
we construct the ratio of two polynomials in n: 

P mi (n) = E" 1 !^ n 1 - 1 ^ 

Qm 2 (n) YT=l d i n ' i ^ 1 ' 

so that P mi (n) /Q m2 (n) match the M—2 ratios a„+i (r) ja n (r) for n = 1, • • • , M — 2. Hence, accounting for 
the arbitrariness of the global normalisation of (|62[) , the complete determination of the two sets of coefficients 
{ci\ i = 1, • • ■ , mi} and {di\ i = 1, • • • , as functions of r implies m\ + rri2 = M — 1. Finally, the auxiliary 
condition on r is obtained by requiring that the last (still unused) ratio au (?") /du—i ( r ) satisfies again the 
n-dependency satisfied by its predecessors, namely that: 

T,?=\* (r) (M-ir 1 = a M (r) 
YTJidi^iM-iy- 1 a M -i(r)- 

The auxiliary condition so obtained is a polynomial in r, the roots of which are candidates to give an 
estimate at order M of r* (noted below r* M ). Notice that, to obtain faster this auxiliary condition, one may 
avoid the calculation of the coefficients q (r) and dt (r) by following the same considerations as those leading 
to with the ADE method. 

At this point, the method potentially reaches the same goal as the ADE and other preceding methods. 
However, according to section 15. 1[ in determining the ratio of polynomials (|62p we have also explicitly 
constructed the function 

F M (x) = r* M ■ mi +iF m2 (ai, ■ ■ ■ ,a mi , 1;/3 1; • • • ,/3 m2 ;a x) , (64) 

in which r* M is the selected estimate of r* , the sets {—ai} and {— are the roots of the two polynomials 
Pm! [ri) and Qm 2 {n) when r = r* M whereas: 

c mi (r M ) . . 

° = d (r* ) ■ (65) 

Now, by construction, Fm (%), has the same truncated series in x as the solution of (|30"|) we are looking 
for. This function is thus a candidate for an approximate representation of this solution. 

It is worth noticing that, contrary to the ADE method, the HFA method does not make an explicit use of 
the conditions at infinity (large x) to determine r*. Only a local information, in the neighbourhood of the 
origin x — 0, is explicitly employed. 

Let us apply the method to the two equations of interest in this paper. 



5.3. Wilson-Polchinski's equation 



5.3.1. Fixed point 

We know that the absolute value of the ratio a n (r*) /a n +i (?**) has a definite value Rwp [given by eq. 
(144[) ] as n — > oo. Consequently, we must consider the ratio (|62j) with mi = m-2 (this implies also that 
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M be odd). In this circumstance, according to section [5.11 the relevant hypergeometric functions have a 
branch cut on the positive real axis (as functions of agx). Consequently the analytic continuation to large 
positive values of x is only possible if o.q < 0. We note also that, according to ([44]) , |ao| should converge 
to 1/Rwp — 0.174774. Finally by considering the large x behaviour directly on (|6Tj) . it is easy to convince 
oneself that the leading power is given by one of the parameters {— c^}, consequently we expect to observe 
a stable convergent value among the ctj's toward the opposite of the leading power at large x of the solution 
looked for. For this reason, instead of the function w (x) of section H.l. li the limit of which is 1 as x — ► oo 
[see (|33)) ]. we have considered the translated function w t (x) = w (x) — 1 which, according to eqs ([5]) and 
([29l |42|) . tends to A*x~ 2 ^ 5 . In this case we thus expect to observe a stable value among the a^s about 0.4 
with the eventual possibility of estimating A*. 

When looking at the roots of the auxiliary condition (|63|) as M varies, we obtain the same kind of 
accumulations about the expected fixed point value r* as shown in figure [1] (with much less points however). 
We can also easily select the right nontrivial solution using the procedure described just above (j39j) . We get 
precisely this excellent estimate with M = 25 and a reduced computing time compare to the ADE method. 
Figure (J6j) shows the accuracies obtained on r* (crosses) compared to the ADE method (open circles). 

Furthermore, the sets of parameters of the successive hypergeometric functions involve two stable quan- 
tities the values of which at M = 25 are: 



a = -0.174 775, (66) 
ai = 0.396 2. (67) 

Those two results are quantitatively and qualitatively very close to the expected values (respectively 
—0.174774 and 0.4 as given just above). 

This clearly shows that the hypergeometric function determined this way provides us with a really cor- 
rect (but approximate) global representation of the fixed point function. This contrasts strongly with the 
numerical integration of the ODE which, due to the presence of the moving singularity, never provides us 
with such an approximate global representation of the solution looked for. 

From (|67|) we have obtained a rough estimate of A* (= 66* /5) by a direct consideration of the value of the 
corresponding function Fm (x) defined in (|64| for some relatively large value of x and we obtain A* ~ —2.6 
what is a sufficiently accurate estimate to serve as a guess in the shooting method. 

We have also tried to determine, using the HEA method, the value A* directly from the "reverse side" 
corresponding to ([43]) . We have not improved the previous biased estimate obtained by ADE (about A* = 
—2.735). We do not understand the significance of this coincidence. We recall, however, that the radius of 
convergence of the Taylor series of u* (y) about y = probably vanishes. This biased result shows again 
that the property of convergence of the Taylor series is crucial for the accuracy of the two methods. 



5.3.2. Eigenvalues 

We have also applied the HFA method to the determination of the eigenvalues. With M = 17, we have 
easily and without ambiguity obtained the following excellent estimates [compare with ([TBI IT8| and tables 
!and[4]: 

^ = 0.649 561774, uj x = 0.655 745 939 3 , w 2 = 3.180 006 53, 
w 3 = 5.912 229 4, w 4 = 8.796 045, ui 5 = 11.8004, 

uj 1 = 1.886 703 839, w 2 = 4.524 390 3 , w 3 = 7.337635. 

These results show a greater efficiency than with the ADE method especially in the determination of the 
subleading eigenvalues. 

It is worth indicating also that, surprisingly enough, we observe again (i.e. as with the ADE method) the 
presence of convergences to the same spurious eigenvalues: 5.8 and 3.8 in the even and odd cases respectively. 
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Fig. 6. Respective approximate number of digits (defined in the caption of figure [2} obtained for r* M with the HFA method 
(crosses) and the ADE method (open circles) for the Wilson-Polchinski fixed point equation. Whereas at a given order M the 
accuracy is similar, a smaller time computing is necessary with the HFA method. 

5.4. Litim's equation 



5.4.1. Fixed point 

Applying the HFA method with the ratio of two successive coefficients a n (k) provides again an accu- 
mulation of roots about the right value of k* given in (|50[) . However, this time, we have encountered some 
difficulties in defining a process of selection of the right root. We obtain the following estimate for M = 21: 

k* ~ 0.409 531, 

which is not bad [compare with ([50)1 ] but not as satisfactory as in the preceding Wilson-Polchinski's case. 

With regard to the transformation (|20[) and the preceding success of the HFA method, it is not amazing 
that the representation of the solution in the Litim case be more complicated than in the Wilson-Polchinski 
case. 

We have already mentioned that, instead of the ratio of two successive terms of the series a n (fc), it is a 
shifted ratio that roughly converges to the finite radius of convergence ([55)1 . As a matter of fact, if we use 
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the ratios 

Qn+3 (fc) 

a n (k) 

instead of the ratio a n+ \/a n without changing the procedurj^l described in section [5?2l then we get a better 
estimate for M = 21 [compare with (|50|) ]: 

k* ~ 0.409 532 737, 

although the convergence properties are not substantially modified. 

Because the case is apparently more complicated than precedently, we do not pursued further the discussion 
of the global representation of the fixed point solution by generalized hypergeometric functions. 

5.4.2. Eigenvalues 

For the eigenvalue problem, a similar difficulty occurs where the right values do not appear as clear 
convergent series of roots. At order M = 17, we get the following estimates: 

v = 0.649 55, ui = 0.657 6, u 2 = 3.20 , u 3 = 5.8 , 
u>i = 1.89 , uj 2 = 4.5 . 

As in the case of the fixed point determination, if instead of applying the method with the ratio of two 
successive terms of the series a n (k) we consider the ratios 

a n +3 (k) 
a n (k) 

then we get better estimates for M = 19: 

^ = 0.649 561774, wi = 0.655 75, w 2 = 3.180 7, uj 3 = 5.905 , 

wi = 1.886 71, #2 = 4.524. 

where the numbers of digits have been limited having regard to the estimated accuracies [compare with ((TSJ), 
table E] (even) and (JTSJ), table H (odd)] . 

6. Summary and conclusions 

We have presented the details of a highly accurate determination of the fixed point and the eigenvalues for 
two equivalent ERGEs in the local potential approximation. First, we have made use of a standard numerical 
(shooting) method to integrate the ODEs concerned. Beyond the test of the equivalence between the two 
equations, already published in [31], the resulting numerics have been used to concretely test the efficiency 
of two new approximate analytic methods for solving two point boundary value problems of ODEs based 
on the expansion about the origin of the solution looked for (field expansion) . 

We have considered explicitly those two methods applied to the study of the two equivalent ODEs. We 
have shown that they yield estimates as accurate as those obtained with the shooting method provided that 
the Taylor series about the origin of the function looked for has a non-zero radius of convergence. 

This is an important new result since, up to now, no such approximate analytical method was known 
to work in the simplest case of the Wilson-Polchinski equation. In the case of the Litim equation the two 
methods converge better than the currently used expansions (usually referred to as I and II in the literature, 
see e.g. [[II]])- Our results support concretely the conclusions of [115] which indicated that the high field 
contributions were important in the Wilson-Polchinski case whereas they were less important in the Litim 
case. 

3 Notice that the procedure does not define some generalized hypergeometric function of x 3 . This would have been obtained 
by considering separately three series in the original series. Then a combination of three generalized hypergeometric functions 
would have represented the solution looked for. 
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The first of the two methods relies upon the construction of an auxiliary differential equation (ADE) 
satisfied by the Taylor series at the origin and to which is imposed the condition of the second boundary (at 
infinity) [15) . 

The second method (HFA) is new. It consists in defining a global representation of the solution of the 
ODE via a generalized hypergeometric function. The HFA method provides the advantage of yielding a 
global (approximate) representation of the solution via an explicit hypergeometric function. 

In both cases it is possible to obtain easily (with few terms in the field expansion) rough estimates of the 
solution which may be used as guesses in a subsequent shooting method. 

The procedures may be applied to several coupled ODEs as shown in [15] for the ADE method. Hence, we 
hope that the present work will make easier and more efficient future explicit (and ambitious) considerations 
of the derivative expansion of exact renormalisation group equations. 
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